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ABSTRACT 



We calculate the "observed at infinity" image and spectrum of the accretion structure in Sgr A*, by modelling it as an optically thin, 
constant angular momentum ion torus in hydrodynamic equilibrium. The physics we consider includes a two-temperature plasma, a 
toroidal magnetic field, as well as radiative cooling by bremsstrahlung, synchrotron and inverse Compton processes. Our relativistic 
model has the virtue of being fully analytic and very simple, depending only on eight tunable parameters: the black hole spin and the 
inclination of the spin axis to our line of sight, the torus angular momentum, the polytropic index, the magnetic to total pressure ratio, 
the central values of density and electron temperature and the ratio of electron to ion temperatures. The observed image and spectrum 
are calculated numerically using the ray-tracing code GYOTO. Our results demonstrate that the ion torus model is able to account for 
the main features of the accretion structure surrounding Sgr A*. 
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1. Introduction 

The Galactic Centre is one of the most interesting regions for 
scientific investigation because it is the closest available galac- 
tic nucleus and, therefore, can be studied with resolutions that 
are impossible to achieve with other galaxies. The radio source 
Sagittarius A* (henceforth Sgr A*) is associated with the super- 
massive black hole at the centre of our Galaxy. Its mass, M = 
4.31 x 10 6 M o , has been determined in particular from the com- 
plete and highly eccentric orbit of the star S2, which passes 
Sgr A* in its pericentre at a mere 17 light-hours, i.e. 1400 r$, 
where rs = 2GM/c 2 is the Schwarzschild (or gravitational) 
radiu s dSchodel et al.l 120021: JBower et all 120041 iGillessen et al] 
120091) . Given the Sgr A* distance of 8.33 kpc, the Schwarzschild 
radius corresponds to an angular size of around 10 //as in the 
sky, making the Galactic Centre black hole an ideal candi- 
date for near- future microarcsecond interferometric technolo- 
gies dPaumard et al] 20081: Eisenhauer et al.l 120081: lEisner et al] 
l2010tlBroderick et al.ll201 lh ~ 



In this context, several authors have calculated images 
of various theoretical models of accretion structures around 
Sgr A*. Numerical simu lations reach fr o m thre e-dimensional 
non-r elativistic MHD dGoldston et alJ 120051 
two-dim ension al 
201 Oh 



Huang et al. 
GRMHD dMoscibrodzka et al. 
and three-dimensional GRMHD 



_ iHilburn et aL ._ 
dDexter et al.ll2010UDibi et al.ll2012h . The size and shape of the 
"black hole silhouette" cast by the black hole on the accretion 
structure, is determined by the photon orbit and depends thus 
only on black hole mass and spin. Therefore, with the known 
black hole mass in the Galactic Centre the spin could be 
estimated from fitting the calculated size and shape of the 
shadow to the size and shape observed. The looks (spatial 
extent, shape, brightness distribution, etc.) of the entire image, 
however, depend not only on black hole mass and spin but also 



on the details of the accretion structure around Sgr A*. Most of 
them are uncertain, e.g., the chemical abundance, the radiative 
processes, the inclination. They have not been sufficiently 
examined yet, and a wide range of the relevant parameter space 
remains unexplored. 

The accretion structure around Sgr A* is most probably 
a radiativel y inefficient a d vectio n dominated accretion flow 
(ADAF), as iNaravan et al.l d 19951) rather convincingly demon- 
strated by means of spectral fitting. In the following, an ap- 
proximate analytic torus model asso ciated to the famil y of 
Polish doughnuts introd uced by lAbramowicz et all d 19781) and 
iJaroszvnski et a l. (1980) is constructed for a radiatively ineffi- 
cient accretion flow to describe the accretion structure in Sgr A*. 
The model depends on a few tunable parameters which can be 
changed and adjusted so that they can cover the whole parame- 
ter space of the problem. The general model of the source, and 
the ray-tracing tailored especially for it, may as well be applied 
to spectral calculations of other black hole candidates, both op- 
tically thin (ion tori, to model e.g. some spectral states in X-ray 
binaries) and optically thick (Polish doughnuts, to model e.g. 
some luminous AGNs and ULXs). 

IShapiro et al.l d 19761) designed a hot, two-temperature disc 
model to describe the strong X-ray emission observed in Cygnus 
X- 1 , which was too hot to be understood in terms of the standard 
model. They found an optically and geometrically thin solution 
branch in which ions and electrons are weakly coupled, having 
different temperatures and being in energy balance. This solu- 
tion, however, is thermally unstable. 



Ion tori were proposed inde pendently b y Ichimarul (Il977l) 
(though not by this name) and iRees et al.l d!982l) to model a 
gas flow that emits little detectable radiation but is at the same 
time able to power radio jets, both in galactic sources and mi- 
croquasars. Ion tori are geometrically thick, gas (ion) pressure 
supported spheroidal structures, located at the inner regions of 
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accretion flows. They are extremely optically thin and very ra- 
diatively inefficient. These ion tori consist of a fully ionised 
plasma, hence their name, in which protons and electrons are 
thermally decoupled (each following its own temperature distri- 
bution) and not in local energy balance. They are assumed to 
be low-a accretion flows threaded with magnetic fields, so that 
the ordinary molecular viscosity is suppressed in favour of mag- 
netic/turbulent viscosity]. Due to their vertically extended shape, 
H/r < 1 , ion tori naturally create a pair of funnels through which 
magnetic flux can collimate and escape. Ion tori only exist in 
the sub-Eddington a ccretion regime, M < M CI \t < ^Edd (see 
e.g. Ilchimarulll977t iRees et all 1 19821) . where radiative cooling 
is inefficient enough that ions stay sufficiently hot (near virial) 
during the whole inflow time, i.e., tag » facer. Otherwise, at 
higher mass accretion rates, cooling causes the puffed up in- 
ner spheroidal region to gradually become opaque and deflate 
to a standard thin disc. As a consequence of fdiff s> f accr , the 
ions carry most of the energy with them into the black hole. 
Instead of "the ions are not cooling down fast enough" one 
may also and equivalently say that the disc is locally "advec- 
tively cooled" because energy is carried away by ions. Therefore, 
this branch of solutions was rename d to "Advection Dominated 
Accretion Flow" (ADAF) (see e.g INaravan & Yilll994t 119951: 



lAbramowicz et aljll995tlEsin et all 19971) 

In this work, we model Sgr A* by io n tori very much similar 
to these described by IRees et al.l d!982l) . They are based on the 
mathematical description introduced by Pac zyriski and collab- 
orators in their studies of Pol ish doughnuts (lAbramowicz et all 
ll978HJaroszvnski et aT1ll980t) . Although Polish doughnuts have 
some of their important characteristics directly opposite to these 
of ion tori (they are optically thick, radiation pressure supported 
and correspond to super-Eddington accretion rates, M » Mbaa) 
both classes of tori have the same equipotential structure and 
very similar dynamical properties. In particular, they both have 
the "Roche lobe", i.e. a critical equipotential that crosses itself 
along the "cusp" at r — r cusp (see Fig. [TJ. Roche lobe overflow 
causes the dynamical mass loss from the torus to the black hole, 
with no need of help from viscosity. Thus, the accretion flow 
at radii r < r cusp is regulated by the black hole strong grav- 
ity and not by viscus processes. The cusp should be regarded 
as the inner edge of the disk; for small accretion rates it cor- 
responds to ISCO, for higher accretion rates the cusp locates 
closer to the black hole. Both ion tori and Polish doughnuts are 
dynam ically unstable with respect to the Papaloizou & Pringle 
(11984 instability. However, as proved bv lBlaesi (1 19871) . this ii> 
stability is suppressed by the Roche lobe overflow (see also 
INaravan & Pophamlll993l) . 

These properties, shared by ion tori and Polish doughnuts, 
are genuine and robust also for much more general toroidal 
structures around black holes. In particular, they do not depend 
on the angular momentum distribution inside the torus. It is cus- 
tomary to assume that angular momentum is constant, I — {q — 
const, as this leads to remarkably simple final analytic formulae. 
This is an assumption adopted here. In follow-up papers we will 
relax this assumption and calculate observed properties of ion 
tori with 6 + con st. Dynamical mo dels for such tori have been 
calculated e.g. bv lOianetall d2009l) . 

The paper is structured as follows. In section [2] we construct 
the hydrodynamical torus and in section [3] and section |4] we dis- 
cuss the radiative properties. Section |5]shows the resulting spec- 
tra in relation to broadband data of Sgr A* and a series of im- 
ages. The presented spectra are neither chi-squared fits nor best 



guesses, they merely illustrate the performance of the ion torus 
model. Section[6]presents the conclusions. 



2. The geometry of fat tori 

2.1. Kerr metric 

In Boyer-Lindquist spherical coordinates (t, r, 9, <p), geometrical 
units c = 1 = G, and signature (-, +, +, +), the Kerr metric line 
element has the form: 



ds 2 = - 1 



2Mr\ , AMra 

sin 2 



- 1 At 2 sin 2 6 At dd> + — dr 2 

a a A 



2Mra sin" #\ , , 
+Zd6 2 + {r 2 + a 2 + \sm 2 6d(f> 2 , (1) 



where M is the black hole mass, a = J/M its reduced angular 
momentum (J being the black hole angular momentum), S = 
r 2 + a 2 cos 2 6 and A = r 2 - 2Mr + a 2 . 

Circular orbits in Kerr spacetime are obtained from the time- 
like geodesies of the metric (Q]i for which the 4-velocity takes 
the form uP - («', 0, 0, u^~). The two constants of geodesic mo- 
tion, associated respectively with stationarity and axisymmetry, 
are then expressible as 



& - -u t - —u (gn + £2 
L - «0 = u' (g tl p + Q 

where 
D.s — 



(2) 
(3) 



(4) 



is the angular velocity with respect to a distant observer. The 
specific angular momentum is defined by 



& Ut 



(5) 



At circular geod esic orbits, the ang ular momentum has its 
"Keplerian" form dBardeen et al.lll972l) : 



Mr,a) = 



M l ' 2 (r 2 -2aM l ' 2 r l ' 2 +a 2 ) 
r 3/2 - 2Mr 1 ' 2 + aM 1 ' 2 ' 



(6) 



Here and in the remainder of the article, we consider M a fixed 
parameter, so that we can take l^ to be a function of (r, a) only. 
Important circular orbits in the Kerr metric are the 
marginally stable and ma rginally bound ones, corresponding to 
the following values of r dBardeen et al.lll972l) : 



r ms ( fl ) = M [3 + zi - ((3 - zi)(3 + zi + 2z 2 )) 1/2 ] , 

(a) = 2M-a + 2M l/2 (M-a) l/2 . 



r m b 



(7) 
(8) 



where zi = 1 +(1- a 2 /M 2 ) 1/3 [(1 +a/M) l/3 + (1 -a/M) 1/3 ] and 

Zi = (3a 2 /M 2 +Z 2 ) 1 / 2 . On these special orbits, the specific angu- 
lar momentum takes the values 



The a-prescription is applied also in this context. 



4, s (a) = ^K(r ms (a), a) and € mb (a) = ^ K ('"mb(fl), a). 



(9) 
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Fig. 1. Meridional cut through a A — 0.3 (black curve) and a A — 0.7 torus (grey curve) around a black hole (the black semi-circle is 
the event horizon). Left panel: Schwarzschild case. Right panel: Kerr case with a - 0.9 M, the long-dashed semi-circle marking the 
ergosphere. 



2.2. Fluid torus of constant specific angular momentum 

We consider a fluid torus of negligible self-gravitation around a 
Kerr black hole. Modelling the fluid as a perfect one, the stress- 
energy tensor is 



[ /JV 



(e + P^Uy + Pg 



fir. 



(10) 



where P is the fluid pressure and e the fluid proper energy 
density. For a purely circular motion {u^ = m'(1,0, 0, Q)), it 
can be shown that the energy-momentum conservation equation, 



V v 7^ = 0, takes the form 



y 

e + P 



-V„ ln(-w,) + 



i -ne' 



(11) 



with ( related to the fluid 4-velocity components by Eq. (0. 
Assuming that the fluid is barotropic, P = P(e), we introduce 
the enthalpy function 



Jo 



dP 
e + P' 



so that th e left-hand side of Eq. (fTTT i becomes simply V^H. 
Following lAbramowicz et all (Il978l) . we consider models for 
which 6 is constant within all the torus: 



= 4 



(13) 



Accordingly, the equation of motion ( fTTT i reduces to V M H 
-Vf! ln(-M,) yielding 



H - W + const, 

where we have introduced the potential 

Ws-ln(-« t ). 



(14) 



(15) 



From the normalisation relation u /J u f ' = -1, we get 



W(r, 0) 



1 



In 



g u + 2Q.g t4 , + Q 2 



(16) 



(g„ + Q.g tll ,) 2 

where D. should be considered as the following function of (r, ff) 
+ togu 



Q(r, ff) 



(17) 



This last relation is easily derived by combining Eqs. (|2), (0 and 

©. 

lAbramowicz et al] (Il978l) proved that the cusp location 
should be between the mariginaly stable and the mariginaly 
bound orbit, which implies that Eq must obey, 



( ms (a) <{q< t mb {a), 



(18) 



with Ans(fl) and t m b(a) given by Eq. (0. Recast in terms of the 
dimensionless parameter 



(12) A = 



io - Cs(fl) 



^mh(fl) - Ans(a) 

the condition ( fT8l is equivalent to 
0<A< 1. 



(19) 



(20) 



Given the radial dependence (0 of the Keplerian specific angu- 
lar momentum, which takes its minimum at r - r ms (a), and the 
constraint ( fT8l . there are two values of r, r m and r c say, for which 
ZyJj, a) - {()■ One has necessarily 



r mb (a) < r in < r ms (a) < r c . 



(21) 



Since for r — r m or r — r c , the actual angular momentum (q 
is equal to the Keplerian one, the gravity and centrifugal forces 
balance each other, implying V^W = 0. From the first integral 
(fT4l . we have then V^H = 0, or equivalently, via (fT2l . V A ,P = 0. 
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A point where the gradient vanishes can be either a singular point 
(self-crossing of an equipotential) or an extremum. In the present 
case, r - r m is the location where one of the equipotentials of 
W, the so-called Roche lobe, self-crosses, a s shown in Fig. [D 
This g ives rise to a cusp at the torus surface (lAbramowicz et all 
119781) . On the other side, r — r c corresponds to the maximum of 
P. The circle r — r c in the equatorial plane is called the centre of 
the torus. 



2.3. Solution for a polytropic equation of state 

To go further, we assume a polytropic equation of state, 

p = Ke l+U '\ (22) 

where K and n are two constants, n being the polytropic index 
and K the polytropic constant. The total energy density e — p+Tl 
is the sum of energy density p and internal energy n, which in 
the non-relativistic limit, II <& p, reduces to e =* p. Equation ( fT2l 
is then readily integrated, yielding 



H = (n+ l)ln(l +Ke xl "). 



The surface of the torus is defined by P — 0. From Eq. (l22l . this 
corresponds to e = 0, and from Eq. (l23l to H — 0. Therefore, we 
may rewrite the first integral of motion ([14} as 



H = W-W S , 



(24) 



where the constant W s is the value of the potential W at the torus 
surface. Denoting by H c and W c the values of H and W at the 
torus centre, Eq. (l24T i implies 

H C = W C -W S . (25) 

Let us introduce the dimensionless potential 

W(r,6)-W s 

co(r, 6) = . (26) 

Wc - W s 

From Eqs. (l24l i and ( f25l l, we have 

H = H c co. (27) 

Substituting Eq. d23l for //, we get 

ln(l+^e 1/ ")=wln(l+^e c 1/ "), 

where e c is the energy density at the torus centre. Solving for e, 
we obtain 



^K^ir-i]" 



(28) 



At this stage, our torus model is determined by five param- 
eters (in addition to the black hole mass M): the Kerr spin pa- 
rameter a, the dimensionless specific angular momentum A, the 
polytropic index «, the polytropic constant K, and the central 
density e c . From A and a, we evaluate Cq via Eq. ( fT9b . The val- 
ues of £o and a fully determine the potential W(r, 9) according to 
formulae ([T6l-(fT7li. Since we are seeking for a Roche lobe filling 
torus, the value W s of the potential W at the torus surface must 
be set to the value at the Roche lobe (cf. Sec. 



W s = W(r m ,n/2). 



Given W s and H c (deduced from e c via (l23Tl). we determine 
W c by (f25t . We know then entirely the dimensionless potential 



u>(r, 9) as given by formula (1261 . We can thus compute the energy 
density everywhere in the torus by formula (128V 

Note that, by construction [cf. Eq. (f26ll. a> is zero at the sur- 
face of the torus and 1 at the centre. Between < a>(r, 0) < 1 the 
toroidal equipotentials go from the biggest possible, the torus 
surface, down to a single line, a circle at r c . Open equipotential 
surfaces have a>{r, 9) < 0. The size of the Roche lobe tori de- 
pends on the spin. The higher the spin the narrower is the angu- 
lar momentum distribution which sets the location of the critical 
radii, and hence the smaller are the tori. For a fixed value of a, 
the Roche torus is the largest, with an infinite outer radius r out , 
when A - 1, i.e., £q - f m b(a). 

3. Adding physics 

In what follows, we set n = 3/2, which corresponds to the adia- 
batic index y=l + l/« = 5/3 of a non-relativistic ideal gas with 
no radiation pressure. This is consistent for the very optically 
thin medium. 



(23) 3.1. Thermodynamic quantities 



In optically thin gas pressure supported ion tori the radiation 
pressure can be neglected. The total pressure P is expressed as 
the sum of the magnetic and gas contributions, P mag and P gas : 

P = Pmag + Pgas (30) 

The magnetic and gas pressures are assumed to be some fixed 
parts of the total pressure (an assumption that is often made in 
analytic models of accretion structures): 

B 2 



mag 



247T 



BP, 



gas 



(1-/3) P. 



(31) 



Here, B is the intensity of magnetic field. The gas is assumed 
to be a two-temperature plasma, with p\, p e , T- u T e being the 
mean molecular weights and temperatures of ions and electrons, 
respectively. The gas pressure is then expressed as the sum of 
the ion and electron contributions: 

Pgd ^ Pi + Pe ^hJll + Il), (32) 

where £b is the Boltzmann constant and m u is the atomic mass 
unit. Let us write the ion and electron temperatures as follows: 



/(w)yU, 



TWu^gas 



: g(w)jUi 



OT u P g 



(33) 



k B e k B e 

where / and g are linear functions of the equipotential function 
oj. 

The condition is that they are equal at the surface: T\ - T e 
when D — 0, and that they are at some decided ratio £ at the 
centre: ^ T; = T e when u> — 1 . That leads to 

Mi 



/(0) = 



/(I) 






= M, 



= M* 



(34) 



(35) 



(29) T, 



j"e + ml; 

and since Eq. d32l is fulfilled when /(<w) + g(cif) = 1, one easily 
finds g(Q) and g(l) which allows to write explicit expressions for 
the two temperatures, 

(l-B)m u P 



[( 1 - oj)M + a>M{] p, 






Aff) 



k B e 
(l-B)m u P 
k B e 



(36) 
(37) 
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Fig. 2. Ion (dashed/dotted) and electron (solid/dot-dashed) tem- 
perature distribution throughout a X = 0.3 torus (black) and a 
A - 0.7 torus (grey) around a Schwarzschild black hole. We as- 
sume a central ion to electron temperature ratio of £ = 0.1. The 
thin black line marks T v j T . 

It is often practical to express equations in terms of the dimen- 
sionless temperatures 



knT e . 



,,2' 



and 6{ - 



k B Ti 



(38) 



Here m e and m\ are the electron and ion masses. 

Figure [2] shows the temperature distribution i n tori. It is dif- 
feren t from the ADAF self-similar solution (e.g.. lNaravan & Yil 
1 19951) which being self-similar must be a monotonic function 
of radius. Here, the non-monotonicity of the temperature curves 
has its roots in the existence of boundaries, the outer at the outer 
torus radius and the inner at the cusp. The ion and electron tem- 
peratures have a maximal separation (given by £) at the torus 
centre and converge towards either end of the torus. Obviously, 
the temperature in our model is zero at these bo undaries. Note, 
that in the torus interior our curves and these from lNaravan & Yil 
(1 19951) are comfortably consistent. Using Eq. (l32t one can now 
derive the relation for the polytropic constant at the torus centre, 



K = 



lcT e 



(l-£)m u e^jUeAt f 



2 / 3 , 



where T ec is the central electron temperature. The total pressure 

is 



P = 



k R T t 



B-fe.c 



c 5/3 



(l-Au e muA^e c 2/3 ' 
and the magnetic field 



B = 



24^ fc B r e , c e 



5/3-1 1/2 



(l-j3)jUe»JuA<f e 



2/3 



(40) 



(41) 



The thermodynamic relations between T, P and e require the 
knowledge of the mean molecular weight /u linked to the ion and 
electron molecular weights by: 



with 
Mi = 



1 


+ 


1 


A'i 


/V 



nim a 



AX + Y 



1 + 3X + Y/2 



and /j. e - 



1 +Y' 



(42) 



(43) 



where X and Y are the hydrogen and helium a bundances, as- 
sumed to be equal to X — 0.75 and Y - 0.25 as in lNaravan & Yil 
d 19951) . On the one hand this implies the effective molecular 
weight of ions to be [i\ - 1.23 and for electrons /Y e = 1.14, 
and on the other defines the electron and ion number densities, 
n e and «;, with respect to the energy density. In the following, 
we will give general formulae in which the factor n\ is replaced 
by the sum over ion species, 



« = Z Z K' 



(44) 



where Z, and rij are respectively charge and number density of 
the j th ion species. 

3.2. Radiative processes 

We have to complete the set of equations by specifying the phys- 
ical processes which will be involved in the radiative cooling. 
We consider a two-temperature plasma cooled by synchrotron 
radiation, inverse Compton process and bremsstrahlung emis- 
sion. A very convenient and general description of such cool- 
ing processes has been presented bv lNaravan & Yil (1 1995b . We 
closely follow their approximations and procedures, and also 
use their formulae (correcting one typo). For the reader's con- 
venience and completeness, we give here a summary of these. 
Electron-po sitron pair creation and anni hilation is neglected but 
as sho wn bv lBiornsson et al] (1 19961) and lKusunose & Mineshiga 
d 1996b this is justified in most cases of interest. 

The following subsections aim at determining the emis- 
sion and absorption coefficients inside the torus due to 
bremsstrahlung, synchrotron radiation and Compton processes. 
The emission coefficient (resp. absorption coefficient), j v (resp. 
a y ), allows to compute the increment (resp. decrement) of spe- 
cific intensity I v when progressing inside the emitting medium 
by a small distance Al: 



(My = jydl, 

dl v - -a v I v dl. 

The cgs unit of j v is ergcirr 3 s~' ster -1 Hz" 
pressed in cirT 1 . 



(45) 



while a v is ex- 



(39) 3.2.1. Bremsstrahlung 



The rate at which energy is lost due to bremsstrahlung, 
fZ. = dEfr/dtdV, includes emission from both ion-electron and 
electr on-electron collisions (ISvenssonll982llStepnev & Guilberj 
I1983I) : 

/br^/d+Zee- (46) 

The ion-electron bremsstrahlung cooling rate is given by 

/ d = « e «cr T a f m e c 3 F ei (6'e) (47) 

where af - 1/137 is the fine structure constant and the function 
F e i(9 e ) has the approximate form in units of [erg s~'crrr 3 ] 



2 ft. 



1/2 



F*{0 e ) = 4[=f) (1 + 1.781 6t 34 ), 



< 1, 



— [ln(1.123 e + 0.48) +1.5], e > 1. 

27T 



(48) 



In the original formula quoted by IStepnev & Guilberfl (119831) 
there is a number 0.42 instead of 0.48 fsee lNaravan & Yilll995l) . 
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For the ele ctron-electron bremsstrahlung cooling rate 
ISvenssonl ([1982) gives the following formula 



/ee = n 2 e rla f m e c 3 F ee (e e ), 



(49) 



where the function F ee (9 e ) is given in units of [ergs ^m 3 ] and 
has the approximate form 



20 



,3/2 



Fee(0 e )= ^(44-3^ 

x(l + 1.10 e + 6£-1.25 e 5/2 ), 6 e <l, 



(50) 



= 240 e [ln(2?70e) + 1.28], 



> 1, 



where r e = e 2 /m e c 2 is the classical radius of electron and the 
Euler number 77 = exp(-yE) = 0.5616 given by the Euler- 
Mascheroni constant, je- And again, we replac e in the original 
formula 5/4 with 1.28 (see lNaravan & Yilll995l) . 

With the bremsstrahlung cooling rate given above one can 
express the bremsstrahlung emission coefficient as 



■ bi- f - l h I hv 1 



4-7rk B T e \ k B T { 



where h is the Planck constant, the 1 /An factor assumes isotropic 
emission in the emitt er's frame and G is the velo city averaged 
Gaunt factor given bv lRvbicki & Lightmanl(ll986h : 



G = I™* 7 ' 



n hv 



1/2 



n \je hv 



3.2.2. Synchrotron cooling 



k B T e 
hv 



k B T e 
hv 



< 1, 



> 1. 



(52) 



The emission coefficient for synchrotron emission b y a relativis- 
tic M axwellian distribution of electrons writes (see iPacholczvkl 
ri970b 



1 e 4mi e v 



M(x M ), 



(53) 



with a factor l/4n again for isotropic emission in the emitter's 
frame and the fitting function, 



M(x M ) 



4.0505a 

1/ 6 



0.40/? 0.53167' 



1/4 



M 2 



exp(-1.8899xi 1 /3 ), 



where 



xm 



2v 
3v o 2 



VQ 



eB 

2mn e c 



(54) 



(55) 



and the par ameters a, B and y are ta bulated for a range of tem- 
peratures in Mahad evan et al.ld 19961) . The fitting formula is valid 
only for e > 1, i.e. T e > 10 8 K, which is satisfied for applica- 
tions to ion tori because the synchrotron emission that dominates 
the central regions of a torus comes, like in other advective flows, 
from relativistic electrons in the tail of the Maxwellian distribu- 
tion. At lower temperatures or in the outer torus regions emi ssion 
is dominated by bremsstrahlung. Below a critical frequency, v c , 
the synchrotron spectrum becomes self-absorbed. As the flow in 



ion tori is very much akin to a spherical flow, this frequency can 
be obtained by equating the synchrotron emission in a sphere of 
some radius, R [cm], to the Ray leigh- Jeans blackbody emission 
from the surface of that sphere and solving for xm which is then 
substituted into 

3 2 
v c = -Vn 0;x M . (56) 

We assume that at low frequencies, v < v c , the absorption is 
locally given by Kirchhoff's law, i.e., the low-frequency part of 
the synchrotron emission behaves like a blackbody. 

3.2.3. Compton cooling 

The soft bremsstrahlung and synchrotron photons in an ion 
torus filled with a thermal distribution of electrons are (inverse) 
Compton scattered to higher energies. Especially in the central 
regions of the flow this can be an important cooling mechanism. 
There is a probability P that a seed photon of some initial energy, 
E m = hv, is in optically thin material scattered to an amplified 
energy E out = J[E m , where 

1 +46> e + 16ftr. 



(51) f= l-exp(-Tes), 



and 



m 



(57) 



Thus, the energy exchange between electrons and photons de- 
pends only on the electron temperature 6 e and the probability 
that a photon will interact with an electron, which is given by the 
electro n sc attering opt i cal de pth r es = J n e crj Al. iDermer et al.l 
(Il99ll) and lEsin et a l. (1996) derived an approximate prescrip- 
tion for the energy enhancement factor due to Compton scatter- 
ing, which is defined as the average energy change of a seed 
photon: 



(58) 



77 = 1 + 771 + 772 — 



where 



9{3K-V) 





1 - 


-VJ( ' 


= 


3-" 


m, 


= 


-1 


-InP/lnJl 



(59) 



The dimensionless energy is given by x — -^7. As the emerging 
photons cannot gain more energy than the electrons they col- 
lide with have, there is an upper limit on the energy, x < 3 6 e . 
Note that in our simple ion torus scenario we apply, without loss 
of generality, the Thomson cross section rather than the Klein- 
Nishina cross section. 

Comptonised emission is 77 - 1 times the seed photon distri- 
bution. The part of the spectrum which can be Comptonised lies 
between the critical synchrotron self-absorption edge, x - x c - 
hv c /m e c 2 , and x - 3 9 e . 

Comptonisation of bremsstrahlung emission is then given by 

v/73+i / „ x^+n) 



;br,C 




1 



773 + 1 



30, 



(60) 



Note, that the corresponding expression in lNaravan & Yil (119951) 
has a typo which we corrected here. 

Comptonisation of synchrotron radiation which is emitted 
mostly at the self-absorption frequency, v c , is given by 



• sy.C 



j^ = j; y lm-m(xJo e r]. 



(61) 



We calculate Comptonisation to second order and assume that 
photons which are up-scattered to e form a Wien tail. 
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3.2.4. Total cooling 

The total emission coefficient is the sum of all radiative contri- 
butions: 



• _ -br -br,C 

JV — Jy T J v 



•sy .sy,C 
■ Jv + Jv • 



(62) 



For a medium in local thermodynamic equilibrium (LTE) at tem- 
perature T, the emission coefficient (l62l and the absorption co- 
efficient are related by means of Kirchhoff 's law: 



(63) 



B V {T) 



where B v is Planck's law of blackbody radiation. For the typi- 
cal temperatures of the ion tori considered in this paper, a v is 
negligible, thus the absorption can be safely ignored. 

Integration of j v over the whole frequency range gives the 
total cooling rate q~ at each point (r, 6) in the torus. Fig.[3]shows 
the different contributions to the total cooling rate from the equa- 
torial plane, 6 - n/2, of a A - 0.3 torus. 



are known analytically for the ion torus, so the integration is 
straightforward. In particular, the absorption coefficient is ev- 
erywhere. 

Photons are launched from the observer's screen, supposed 
to be spatially at rest at r - 100M, in a given specified solid an- 
gle. The null geodesies are integrated until the torus is reached 
(or until the photon escapes too far from the torus) and the spe- 
cific intensity I v is updated at each step inside the torus by using 
Eq. d45b . A map of specific intensity is thus obtained for a given 
set of directions on sky. It is then straightforward to compute the 
observed specific flux F v by summing these values of specific 
intensities over all directions of incidence. 

Two kinds of G YOTO-computed quantities are presented be- 
low. The spectra (Figs. |4]-|9]l show the quantity v F v for different 
values of frequencies. The images (Fig.fTTJb are maps of specific 
intensity I v . 

Let us stress that the whole C++ code for the ion torus is 
now included in the open source GYOTO code, available at the 
above mentioned URL. 
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Fig. 3. Radiative cooling rates of a A - 0.7 torus around a spin 
a* = 0.9 black hole. The contributions from bremsstrahlung 
(black, dashed), synchrotron (black, solid), Comptonised 
bremsstrahlung (grey, dot-dashed) and Comptonised syn- 
chrotron emission (grey, dotted) are shown as a function of ra- 
dius. 

Radiation that originates close to a black hole is influenced 
by various relativistic effects, such as the gravitational bending 
of light rays, gravitational redshift and Doppler beaming. We 
ray-trace the emission derived in Eq. (l62t from each point (r, 6) 
inside the Roche lobe equipotential back to the observer explic- 
itly including all relativistic effects in light propagation. This is 
described in the next two sections. 



4. Ray-tracing 

We use the General relativitY Orbit Tracer of the O bservatory of 
Parish (GYOTO), presented in lVincent et all d20TTb . GYOTO al- 
lows integrating null and timelike geodesies in any analytical or 
numerical metric. It is optimised to compute geodesies in Kerr 
spacetime efficiently, and is able to integrate the radiative trans- 
fer equation in optically thin media. 

The required inputs for GYOTO are the values of the 
emission coefficient, j v , and absorption coefficient, a v , at any 
point inside the torus. These quantities, defined in section 13.21 



5. The spectra and images of ion tori 

This section presents the spectra and images of an ion torus sur- 
rounding a Kerr black hole. Table [T] gives the reference values 
of the parameters used for the computations. The value chosen 
for the c entral density e c i s standard for Sgr A* (see e.g. the 
values in lYuan et alj|2003l iLiu et al.ll2004t iMoscibrodzka et al.l 
12009b Dodds-Eden et al.ll2010l) . Various parameters are varied in 
the computations: the spin, a, the dimensionless angular momen- 
tum, A, the magnetic to total pressure ratio, B, the electron to ion 
temperature ratio, £, the ratio between central electron tempera- 
ture and virial temperature, and the inclination, i. 

The following sections F5TTI and 15721 show the impact of these 
parameters on the spectrum and image of the ion torus. 



parameter 


value 


a 


0.5M 


,1 


0.3 


P 


0.1 


n 


3/2 


e c [gem -3 ] 


io- 17 


f 


0.1 


T 

1 e,c 


0.02 r vir 


*[°] 


80 



Freely available at the URL http : //gyoto . obspm.fr 



Table 1. Parameters for our reference ion torus. Parameter i is 
the inclination angle of the black hole rotational axis towards 
the line of sight. 7\,i r is the virial temperature. 



5.1. Torus spectra 

Figures [4] to [9] show the impact of various parameters on the ob- 
served spectrum of an ion torus. All parameters that have not ex- 
plicitly specified values are set according to Tab.Q] These figures 
show that the ion torus model is able to account for the general 
features of the observed data, by appropriately tuning the various 
parameters. This includes the X-ray flare "bow tie" which in the 
ion torus model may originate from soft photons that are inverse 



O. Straub et al.; Imaging ion tori around Sgr A* 



Compton scattered by the same population of hot electrons that 
is responsible for the synchrotron emission. Only the flattening 
of the spectrum at low frequencies is never matched, which is 
due to the absence of a non - thermal electron distribution in our 
model, as demonstrated by lYuan et al.l (120031) . A more evolved 
model taking into account this effect will thus be developed in 
future work. 

Figure |4] shows that different spins significantly alter the 
emission. This is due to the fact, that at high black hole spins 
the accretion torus shrinks and moves very close to the black 
hole. Therefore, the flux is being shifted to slightly higher 
energies. Increasing spin enhances synchrotron flux, while 
bremsstrahlung is being softened. 

Figure [5] shows that increasing the central temperature im- 
plies a displacement of the whole spectrum towards higher 
fluxes. The same effect is obtained by increasing the dimension- 
less angular momentum A (Fig. [6]). This is due to the fact that the 
torus puffs up like a balloon when A increases (see Fig.[T]i, thus 
providing more flux. 

Figures [7] and [8] show that increasing the pressure ratio /3 or 
temperature ratio £, has an opposite impact: a higher /? gives a 
higher synchrotron flux while a higher f gives a smaller flux. 
This is due to the fact that increasing /3 increases the magnetic 
field strength, while increasing £ decreases B (see Eq.l4TTi. 

Figure [9] shows that the inclination has little effect on the 
spectrum. Although the number of illuminated pixels is smaller 
at higher inclination, the beaming effect is more important. 
These two phenomena have opposite effect, and the resulting 
spectrum is not changed much. 

Future work will be devoted to studying the possibility of 
constraining the various parameters, in particular Sgr A*'s spin, 
by fitting such spectra to observed data. 
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Fig. 4. Impact of spin on an ion torus spectrum: a — (solid 
blue), 0.5M (dashed red) or 0.9M (dotted green). All other 
parameters are set to their reference valu es listed in Tab, [p 
The black o bserved data are taken from IZhao et all (1200 ll) . 
Zvlka et al. 



(1 19951). 



(120081) for ra d io and sub- 



. , Marrone et al. , 
mm da t a.lTelesco et all(]1996j) . ICotera et al.ldl999hjEckart et ail 
(2006), Schodel et al. (2007) for far- and mi d-infrared data, 
Gen zel et al.l (120031) for near-infrared data and iBaganoff et al.l 
(120031) for the X-ray bow tie. The down pointing arrows refer 
to upper values. Note that we do not present fitted spectra here. 



5.2. Torus images 

Each image is a superimposition of a first-order image (the 
thick distorted ring) and higher order images (the circles cen- 
tred around the black hole) which result from photons that swirl 
around the black hole before reaching the observer. There is 
also a very fine circle of light in each image. It consists of 
photons th at come from a loca tion just outside the photon or- 
bit (see e.g. lBardeen et al.lll972l for a definition). Photons which 
escape from the region inside the photon orbit are so severely 
red-shifted that they create an area of reduced intensity on the 
observer's screen: this is the so- called "black hole silhouette" 
(also known as black hole shadow JFalcke et alfeOOOh . 

The interest of such images is the following. Knowing the 
particular extent and shape of the accreting region from obser- 
vations allows to constrain various parameters of the flow. In par- 
ticular the location and shape of the photon orbit contains infor- 
mation on the geometry of spacetime. Increasing spin displaces 
the thin photon ring in the image plane off centre and a, > 0.9 
introduces in addition distortion. Black hole spin, torus dimen- 
sionless angular momentum and inclination have each a huge 
and characteristic impact on the observed image of the torus as 
depicted in Fig. [TO] This is also revealed in the spectra Fig.|4]and 
[6] The inclination angle, however, does not influence the spectra 
in a commensurate way (Fig.|9]l. 

The required technique of measurement of the black hole 
silhouette is already available thanks to recent progress in mil- 
limetre Very L ong Baseline Interferometry (mm-VLBI) (see 
iDoeleman et alJl2008h . With the near-future sub-mm VLBI in- 
struments direct observation of processes in the range between 
5-40 Schwarzschild radii will become feasible and with it also 




Fig. 5. Impact of central temperature: T /T viT = 0.01 (solid 
blue), 0.015 (dashed red) or 0.02 (dotted green). 



the potential to image radiation from an accretion structure and, 
possibly, to deduce the existence of an event horizon. 
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6. Conclusion and perspectives 

We have calculated electromagnetic spectra and images of ac- 
cretion structures around the central black hole in Sgr A* using 
a simple analytic model of "ion tori " very sim i lar to these de- 
scribed in the well-known paper by iRees et aD (Il982l) . Results 
depend on observationally unknown tunable parameters of the 
model, in particular the black hole spin. The hope is that fit- 
ting these analytic models to observations could practically re- 
strict the allowed parameter range, so that the future sophisti- 
cated MHD simu lations of Sgr A * (similar to recent and most 
advanced ones bv lDibi et al.ll2.012h will be more targeted. 



One interesting question here is how accurate are the simple 
analytic models in comparison with advanced MHD simulations. 
We are working on a comparison of o ur spectra and images with 
these calculated bv lDibi et al.l d2012l) . If the preliminary impres- 
sion that these images and spectra are indeed very similar would 
survive the test, then it would be much easier to fit theory to ob- 
servations, and in particular to eventually measure the black hole 
spin in Sgr A*. 

We plan to expand our models of ion tori to include a more 
complete coverage of the parameter space in conn ection with 
a "ge neral" angular momentum distribution (as in lOian et all 
120091) and to include, directly and analytically, large scale mag- 
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netic fields (as in iKomissar ov 2006). We also plan to cal- 
culate models of ion tori in a non-Kerr background (as ad- 
vocated by Psaltis and co llaborators, e.g., iPsal tis et al. 2008; 
iJohannsen & Psaltisl 120 1 Oh to possibly constrain gravity in the 
strong-field regime of gravity theories alternative to Einstein's 
general relativity. In the context of flares it m ay also be interest- 
ing to examine non-stationary ion tori (as in lAbramowicz et all 
[1981 . 
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Fig. 10. Upper left: Image of the reference ion torus corresponding to the parameters listed in Tab. \T\ as observed by an observer on 
Earth, in inverted colours (the darker, the more luminous). The display is intentionally saturated for a better rendering of the less 
luminous parts of the image. Upper right: same as in the upper left but with the spin parameter increased to a - Q.9M. Lower left: 
same as in the upper left but with dimensionless angular momentum increased to A — 0.6. Lower right: same as in the upper left but 
with the inclination angle decreased to i - 40°. 
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